Dissipation and coherent effects in narrow superconducting channels 
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We apply the time dependent Ginzburg-Landau equations (TDGL) to study small ac currents of 
frequency cu in superconducting channels narrow on the scale of London penetration depth. We show 
that TDGL have t-dependent and spatially uniform solutions that describe the order parameter with 
an oscillating part of the double frequency coexisting with an ac electric field. We evaluate the Ohmic 
losses (related neither to the flux flow nor to the phase slips) and show that the resistivity reduction 
on cooling through the critical temperature Tc should behave as (Tc — T)^ /lo'^. If the channel is 
cut out of an anisotropic material in a direction other than the principal axes, the transverse phase 
difference and the Josephson voltage between the channel sides are generated. 

PACS numbers: 74.20.-z,74.20.De,74.50.+r 



I. INTRODUCTION 

It is a common knowledge that superconductors dissi- 
pate in the presence of the flux flow or, for large driving 
current densities, due to phase slips. It is also known 
that even a small ac current in zero applied field causes 
dissipation when none of the above sources of dissipation 
are present. E.g., the resistive transition to the supercon- 
ducting state recorded with small ac currents has always 
a finite width which for small enough currents and in zero 
field not always can be relegated to the flux flow, sam- 
ple inhomogeneities, or thermal fluctuations. A qualita- 
tive explanation of this dissipation employs the two-fluid 
model with the normal and superfluid densities, n„ and 
Us, constant in space and time.^ 

Ohmic losses in superconductors are absent for small 
dc currents. As was originally argued by Landau for su- 
perfluids, the flow of quasiparticles is stopped by the lat- 
tice (phonons) or by impurities and does not contribute 
to the current, whereas the creation of new excitations 
is prohibited by the gap in the quasiparticle spectrum. 
The situation is different for ac currents. During the ac 
period 27r/a;, the normal part of the Fermi- liquid does 
not stop completely and, therefore, causes Ohmic losses. 
When CO significantly exceeds the phase relaxation rate 
tJ~^, but still is small relative to the normal carriers re- 
laxation rate , the Ohmic losses should approach their 
normal limit ja. 

These results were obtained within microscopic the- 
ory, see e.g., Ref.|2- In this article we show that for 
low frequencies the TDGL offers a general and simple 
method to approach the dissipation problem near the 
transition point without specific assumptions on the dis- 
sipation mechanism. We show that if superconducting 
wires (channels) are thin compared to the London depth 
A and the ac current can be taken as uniform, the order 
parameter acquires a part oscillating in time with the fre- 



quency Ilo where to is the current frequency. The order 
parameter modulus stays constant in space since no vor- 
tices or phase slips are assumed to exist in zero applied 
field and for sufficiently small currents. One can say that 
there is a periodic exchange between the superfluid con- 
densate and the normal excitations, accompanied by an 
ac electric field E. In general, the phase shift between 
the field and the current depends on relative values of 
WTA and LOT J with ta and tj being the relaxation times 
for the order parameter and for the phase (i.e., for the 
current). As a consequence, the dissipation depends on 
these parameters, too. 

In anisotropic superconducting channels, the ac cur- 
rents flowing in any but the principal crystal directions 
cause the electric field to have a component perpendicu- 
lar to the current, i.e., across the channel. This is due to 
the anisotropy of the superconductor in use and due to 
anisotropy of the normal conductivity. We show that for 
LOT J <C 1 the transverse field is caused by the inherent for 
anisotropic superconductors transverse phase difference. 
This offers a relatively simple probe of existence of this 
phase difference which has been recently predicted^ 



II. ISOTROPIC CASE 



To set notations, we start with the first GL equation: 



-^^n^i/; = V (1 



l^lV^o) ■ (1) 

Here, ^ is the coherence length, 11 = V+27riA/(^o with A 
and ^0 being the vector potential and the fiux quantum. 
For the order parameter written as ip = f e^^ , we have 
Hip = e^^{Wf + iPf) where P is proportional to the 
gauge invariant vector potential 

Q = 0oVx/2^ + A-(/.oP/27r. (2) 

Equation ^ contains e*-^ on both sides. After can- 
celling this factor and separating real and imaginary 



2 



parts, one obtains for the real part: 

-e'(vV-/p') = /(i-/V/o). 



(3) 



The imaginary part coincides with divji = div/^P = 0. 
The gauge invariant form of TDGL involves the scalar 
potential (pii^ 



d_ 

dt 



^^^)^ = V^(^i-^Ke'nV, (4) 



where t\ is the order parameter relaxation time. Sepa- 
rating real and imaginary parts we have: 



dt 



e{y^f-fp'), (5) 

-TAc/$ - V/-Q + V-(/g), (6) 
where $ = — (0o/27rc)9tX- 



III. SPATIALLY UNIFORM SOLUTIONS OF 
TDGL 

We are interested in coordinate independent solutions 
f{t) and Qi{t). The system (jsj, ^ then takes the form: 



TA du 



u) - fP^ u , 



0. 



72' (7) 



(8) 



This should be complemented with equations for the cur- 
rent. A uniform current J consists, in general, of normal 
and superconducting parts: 



J ^aE 



2e2 , 



(9) 



where E is the electric field directed along the channel 
and a is the conductivity for the quasiparticles flow. We 
aim to describe the system response to ac currents; then a 
is in general w-dependent. If however the frequencies are 
bound by inequality a;T„ <C 1 with r„ being the scattering 
time for the normal excitations, one can consider cr as a 
real w-independent quantity. 

The electric field is expressed in terms of gauge invari- 
ant potentials: 



E ■ 



cdt 



dQ 

cdt ' 



so that the total current is 



IfdQ 2e2 



(10) 



(11) 



At a given current, Eqs. Q and Ijllll form a complete 
system for two functions u{t) and Q{t). It is convenient 
to introduce dimensionless vector potential 



27rC 



(12) 



and to measure the current density in units of the de- 
pairing value:^' j = J / Jgl with 



Jgl — 



C0O 



Also, we use the so-called "current relaxation time' 



Ma 
2^ 



n.. 



(13) 



(14) 



In our view, a better term for this quantity is the "phase 
relaxation time" which we use in what follows, but we 
keep the standard notation tj. When T ^ Tc, tj oc 
1/iTc- T) and so does ta.^ 

Then, the system of equations to solve takes the form: 



-Tj q-uq = j, 



(15) 
(16) 



where dots stand for d/dt. 

We are interested in calculating the system response 
to ac currents J = Jq cosujt with amplitudes Jo Jgl, 
i.e., for j -C 1. In this situation the order parameter u 
is close to unity and q <ti I. The system to solve can be 
simplified (u « 1 — w, v ^ 1): 



TA v/2 + v = q^ , 

Tj <i + Q — ~ jo cos LUt . 



(17) 
(18) 



The second equation here is linear, moreover, it is decou- 
pled from the equation for v and is easily solved. The 
solution consists of a transient part depending on ini- 
tial conditions (the general solution of the homogeneous 
equation) and the long time asymptotics of our interest 
(the particular solution of the full equation). The lat- 
ter can be readily found by looking for q of the form 
A sin Lut + B cos cut: 



q(t co) 



jo{ujTj sinut + cos ivt) 



1 



(19) 



This can also be written in a more familiar form: 



Jo 



sm{ujt + a) , tana = . (20) 

LUTJ 



In the following we are interested only in the station- 
ary long time asymptotics and omit the subscript oo. 
Substituting the obtained q in Eq. IjlTI) we can find the 
long time asymptotics for v. To this end, we look for 
V = Vq + vi COS 2wi -|- V2 sin 2Lut and obtain: 



vo 



Vl 



V2 



2(l+^Vj) 



jg(i-cjV3-2cjVjTA) 

2(1 + WV2)2(1+ 0,2^2) 
2(1+^2^2)2(1+^2^2) 



(21) 
(22) 
(23) 
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This yields: 



j1 



j1 



2(1- 
tan P ~ — . 

V2 



sin(2wi + /?) , (24) 
(25) 



Thus, in the stationary state reached when t ^ 
niax(rA, Tj), the order parameter has a part osciUating 
with frequency 2uj near the average value given in the 
first two terms of Eq. (|24|l . Clearly, the frequency dou- 
bling is due to the order parameter independence on the 
current direction. The zero frequency limit of Eq. I|24|) 
coincides with the known GL result for the order param- 
eter suppression by a dc current: u = 1 — Jq. 

Oscillations of the order parameter per se are difficult 
to measure. This is not the case for the electric field E 
and the dissipation density W = JE. The field E of 
Eq. H1U|I in the stationary long time state is 



E = 



<Poq ^ _0o 



cos(wt + a) . (26) 



The dissipation averaged over the oscillations period fol- 
lows: 



Tj 



HI 



(27) 



It is worth noting that for small currents both the electric 
field and dissipation are not affected by the order param- 
eter relaxation time ta- For ujtj <C 1, the field E on lu 
and W (X Lu'^; they are w-independent for ujtj ^ 1. 

Since tj diverges when T ^ Tc, we obtain in this limit 
the dissipation in the normal state W = Jg/2(T, as ex- 
pected. Expanding W of Eq. (|27|l in the small parameter 
1/lu'^Tj and keeping the first correction we obtain 



W : 



2(7 



1 - 



(28) 



While looking at the T-dependence of the dissipation 
near Tc, it should be noted that the quasiparticle conduc- 
tivity a = Une^Tn/M decreases linearly in (Tc — T) due 
to a decrease of the normal density Un- This causes an 
initial increase of W, which can be considered as mani- 
festation of well-studied coherence effects in electromag- 
netic absorption. However, for frequencies of our interest 
<^Tn <C 1, the "bump" (the maximum) in the dissipation 



AW ~ ( Jq/ct) w2r2 is situated at T « Tc(l 



T-n), i-e., 



very close to Tc- Out of this narrow temperature domain 
the dissipation reduction on cooling through Tc should 
behave as {Tc~T)^/uj^. 



IV. ANISOTROPIC CHANNEL 

In isotropic superconductors in the presence of per- 
sistent currents, the gauge invariant gradient VO = 



Vx -I- 2TTA/(j)o is directed along the current. Recently 
Kogan and Pokrovsky showed that in anisotropic super- 
conductors the transverse phase difference may appear if 
the driving current does not point in any of the principal 
crystal directions."^ In particular, this situation is real- 
ized in current carrying channels cut out of anisotropic 
crystals with a long side in any but a principal direction 
and which are narrow on the scale A. One of the possible 
ways to observe the transverse phase is to measure the 
voltage V generated by time-dependent phase difference 
according to the Josephson formula V — {h/2e)dAQ/dt. 
This can be achieved by driving an ac current through the 
said channel, a simpler possibility to observe the trans- 
verse phase than that suggested in Ref. [sl 

In the static case, the supercurrent density is given by 
Ji = 2ehM^^\A\'^Pk, where Mij is the superconducting 
mass tensor; the summation is implied over repeated in- 
dices. It is convenient to introduce the normalized inverse 
mass tensor = M~^M with M = {MaMbMc)'^/^] 
then the eigenvalues are related by ^a^^b^^c = 1- In 
the uniaxial case, /l(^/Xc = 1, the inverse masses can 
be expressed in terms of a single anisotropy parameter 

72 = ^ia/|Jl,c■■ IJ-a = 7^/^ Mc = 7""^^- 

In the coordinates of Fig^ the components are: 



= 7"*/'(7'sin2 + cos2 0), 

fixy = 7"*^^(l-7')sin0cos0, 
whereas fXzz = fJ-b = 1^^'^ and ji^x = l^zy = 0. 



(29) 




FIG. 1: The superconducting channel with a long dimension 
along X, the direction of an ac current. The crystal axis b is 
directed along z. The other crystal axes are in the xy plane 
with a misalignment angle 6 between c and x. 

To describe i-dependent situations, we again employ 
the TDGL model, the generalization of which to the 
anisotropic situation is straightforward: one has to re- 
place the operator ^2n2 in Eq. Q with fiikHiHk (see, 
e.g., Ref.^. Then, we employ the same procedure as in 
the isotropic case to make the model dimensionless. The 
scalar quantities ^ and A now have meaning of averages 
(Cq6Cc)^/^ and (AaAbAc)^/^, respectively. 

As a result, the system (|15|l . (|16|l is replaced with 

TAu/2 = u-u^ ~ tJ.ikqiqk u , (30) 
TjSikqk + unikqk = -joSix cosujt . (31) 

The tensor Sik — Oikl'^ is the normalized conductiv- 
ity with o = [paOhO^^I'^ ■ As with the mass tensor, 
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we can introduce for the uniaxial case the conductivity 
anisotropy parameter 7^ = aa/(Jc so that Sa = Sf, = 
and Sc = ^a'^^'^ ■ With these definitions, the components 
of Sifc are given by formulas (|29|l in which 7 is replaced 
with 7(j. 

For small currents, we have for v — \ — u: 

-T/^v/2 = - ^.ikqiqk , (32) 

TjSikqk + y'ikqk = ~jo5ix cos ut , (33) 

As in the isotropic case, the equation for qi is decoupled 
from that for v. One can look for a particular solution of 
Eq. H33|l in the form qi = Ai sin ujt + Bi cos ujt to obtain 
a linear system of equations for Ai, Bi. 

Perhaps, the easiest is to deal with Eq. Ip^ in the crys- 
tal frame (a, b, c) where all material tensors are diagonal. 
In this frame, the equation to solve reads: 



TjSaqa + fJ-aqa = -joa COS Ujt , a^a,C; 

joa = - jo sin 6 , joc = jo cos 6 
The long time asymptotics is easily obtained: 

joaSmiujt + Pa) , „ ^J-a 



qa 



tan Pa 



LUTjSa 



(34) 



(35) 



The electric field components in the channel frame (x, y) 
and the dissipation read: 



In the linear approximation in the small wtj, the time 
averaged dissipation is absent; the energy during each 
period is pumped from the condensate to the normal ex- 
citations and back in equal amounts. The electric fields 
are: 



Ey — Ex 



(39) 



The conductivity tensor Sik does not enter these expres- 
sions. One may say that these electric fields are due to 
the t-dependence of the phase differences (the Josephson 
relation mentioned above). In particular, the very fact 
that the transverse field Ey is a proof of existence of 
the transverse phase. Hence, measuring the transverse 
and longitudinal voltages on a channel similar to the 
shown in Fig.^ one can, in principle, verify Eq. (|39f) and 
therefore observe the transverse phase difference. 

3. LOT J ^ 1, the situation taking place in particular 
when T Tc- The dissipation of Eq. H36|l reduces to the 
form similar to that of the isotropic case: 



1 - im 



(40) 



with 



Ex — 



Ey = 



W = 



sm 



' cos((jji 4- Pa) cos^ 9 cos{ujt + Pc) 



t4 



4>ojo^ sin 29 
47r^c 



cos{u!t + Pc 



COs{iLjt + Pa) 



Sn sm 



Sr cos 



(36) 



Clearly, these expressions have the correct isotropic limit. 
It is instructive to consider a few limiting situations. 

1. For a dc current (a; — 0), Eq. H32|l gives v = ^likqiqk 
whereas Eq. yields qi = —jofJ-xl- We then have 



— 1 jol^xx I 



(37) 



so that the order parameter suppression by a dc current 
depends on the current direction. We will not write 
down a cumbersome expression for u in the general case, 
but the physics here is the same as for the isotropic case: 
the order parameter has a small part oscillating with 
frequency 2w. 

2. u!Tj ^ 1. This situation corresponds to temper- 
atures not too close to the critical temperature because 
Tj —> 00 when T — > Tc. We have: 



— J( 
2a 



W^'^u^t] f ^sin^e+^cos^' 



(38) 



sm p -I- 7^ cos^ ( 



cos^ 9 



(41) 



Ki W — J'^a^xl'^ is the normal state dissipation. 
Thus, on cooling through Tc, the resistivity drop should 
behave as (Tc — T)"^ juP' with an angular dependent coef- 
ficient. It should be noted that 70- may exceed substan- 
tially the superconducting anisotropy 7 causing a strong 
angular dependence of r\. 

We have also performed the linear stability analysis of 
our solutions of TDGL to show that the homogeneous 
solution is stable unless the current reaches the magni- 
tude of the order of Jql- One can argue that vortices 
might be generated near the boundaries at smaller 
currents thus destroying the uniform time-dependent 
states. Without going to a detail discussion of this 
restriction, we note that our solutions for small currents 
are certainly stable. 
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